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We discuss the use of super-fermion formalism to represent and solve quantum kinetic equations 
for the electron transport problem. Starting with the Lindblad master equation for the molecule 
connected to two metal electrodes, we convert the problem of finding the nonequilibrium steady state 
to the many-body problem with non-Hermitian Liouvillian in super-Fock space. We transform the 
Liouvillian to the normal ordered form, introduce nonequilibrium quasiparticles by a set of canon- 
G ' ical nonunitary transformations and develop general many-body theory for the electron transport 

through the interacting region. The approach is applied to the electron transport through a single 
level. We consider a minimal basis hydrogen atom attached to two metal leads in Coulomb blockade 
r~^ ' regime (out of equilibrium Anderson model) within the nonequilibrium Hartree-Fock approximation 

as an example of the system with electron interaction. Our approach agrees with exact results given 
by the Landauer theory for the considered models. 
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I. INTRODUCTION 



Recently, there have been a significant progress towards theoretical and experimental understanding of electron 

transport through single molecules ii"— The research has lead to discovery of interesting physical transport phenomena 

j^ . in molecules such as Franck-Condon blockade,^ nonequilibrium Kondo effect,™"™ negative differential resistance;^ 

vibronic effects and local heating,^*' as vifell as svi^itching and hysteresis. ^^ Being a nonequilibrium system of strongly 

r ■ correlated electrons molecular junstions present many interesting and challenging physical problems. Because of the 

■ O ■ very large current density (for example, if we consider a typical experimental current of a few fj,A across the molecular 

C , junction of one atom wide, then the current density would be much larger than in usual mesoscopic devices)** and 

high inhomogeneity of electron density (such as, for example, cusps on nuclear positions or electron concentrations 

along chemical bonds), the electron-electron correlations play a pivotal role in electron transport through molecules. 

Much of the theoretical and computational studies of transport properties of molecular nanostructures are based 

f^ ' on Landauer theory which defines the conductance from the electron transmission probability.*^ The transmission 

probability is usually computed by nonequilibrium Green's functions (NEGF).*'^ Most of the currently employed 

first principles electron transport calculations combine NEGF or scattering theory type approaches with ground 

\^ state density functional theory (DFT)J^"— There are many diverse practical implementations of NEGF-DFT and 

T^ ' it has been applied to various molecular j unctions J^^—iiSii^"— The achieving of the agreement between these first 

principles electron transport calculations and experiments has been elusive target for the last decade: the theoretically 

J-^ ■ predicted current is systematically orders of magnitude too large i^^i*^'^^" — Since reliable and accurate experimental 

f~^ measurements of single molecule transport properties are now available ;2i2£ we cannot anymore blame poorly controlled 

experiments for these discrepancies. All practical applications are limited to mean field, and, in spite of some recent 

efforts^"— electron-electron interaction effects are difficult to include into practical calculations. Therefore, we think, 

it is very important at the current stage of research to move beyond currently employed NEGF-DFT scheme towards 

p\. • successful highly accurate methods for quantum chemical electronic structure calculations, such as, for example, 

H coupled cluster theory^^'^^ and configurational interactioui^^ To succeed we need to develop a formalism which enables 

■ - - us to extend ideas of these sophisticated methods for equilibrium electronic structure calculations to nonequilibrium 

systems, namely, we ideally would like to reduce the problem of electron transport calculations to the eigenvalue 

problem with large, non-hermitian, but nevertheless finite matrix. This paper presents our first step in this direction. 

The quantum kinetic equations provide a more suitable language for molecular type electronic structure calculations 

than NEGF or scattering theory based approaches. Recent progress in the development of kinetic equations for the 

electron transport presents a significant promises in this direction)^"— Although to fully engage transport kinetic 

equations with the high level electronic structure calculations we need to transform them to the familiar language 

of second quantization, creation and annihilation operators, normal ordering, Fock space and vacuum. In this paper 

paper we develop a systematic scheme to convert and solve quantum kinetic equations in the language of the advanced 

quantum chemistry. Our derivation is based on the Lindblad master equation,*- but the approach can be readily 

extended to other types of kinetic equations. 

We would like to distance the formalism presented in this paper from the recent work of one of the authors. ^^ There 



was proposed a theoretical approach to the electron transport problem based on methods of equilibrium thermofield 
dynamics. The approach used the square root of the density matrix and employs unitary canonical transformations 
to build nonequilibrium Fock space. Because of that the main restriction of the method is the necessity to use 
the full space Liouvillian, which includes the bath continuum states, throughout the calculations. Contrary, in the 
current paper we develop a set of nonunitary canonical transformations which enables us to lift the above mentioned 
restriction. 

The important remark on the notation used in the paper is the following: only creation/annihilation operators 
written with letter a (such as for example aka and cli.^) are related to each other by the hermitian conjugation; all 
other creation (c^, 5^, S etc) and annihilation (c, &, d etc) operators are "canonically conjugated" to each other, i.e. 
for example h^ does not mean {b)'^ although {6,6^} — 1. 

The rest of the paper is organized as follows. In Section II, the aspects of the super-operator formalism relevant 
to the electron transport problem are described. We also demonstrate in this section how to solve the Lindblad 
master equation in the super-fermion space for a single impurity connected to one thermal bath. Section III presents 
the general formalism to work with the Lindblad master equation in the super-fermion Fock space for the electron 
transport through the interacting region. In section IV, we apply the formalism to the electron transport through 
single level and through molecule in Coulomb blockade regime. In this section we also compare the results with the 
Landauer theory. Conclusions are given in Section V. We use natural units throughout the paper: h = ks ^ \g\ = 1, 
where — |e| is the electron charge. 

II. LIOUVILLE SPACE AND SUPEROPERATOR FORMALISM 

In this section we describe the algebraic structure of the super-fermion space relevant to the electron transport 
calculations. We also establish the notation and terminology which are used throughout the paper. This section is 
partly based in the formalism of superoperators developed by Schmutz in the context of real-time Green's functions 
method^. Some aspects of the superoperator nonequilibrium Green's function theory with application to electron 
transport have been recently discussed by Harbola and Mukamel.^^ In contrast to Harbola and Mukamel's work, the 
present approach does not aim to re-formulate Keldysh NEGF in terms of super-operators and uses the formalism 
to solve quantum kinetic equations for the electron transport problem. We also use some ideas laid down by Prosen 
in his works on "third quantization" techniquC) ^^'^^ where the similar method was used to solve master equation for 
non-interacting open Fermi systems with quadratic Hamiltonians and one dimensional spin chains. 

A nonequilibrium quantum system can be described in terms of density matrix p(t) which satisfies the Liouville 
equation: 

zm = [H,p{t)]. (1) 

The average value of an operator at particular moment t is given by 

{A{t))^Tr{p{t)A). (2) 

For many-particle quantum systems the density matrix p{t) and the Hamiltonian are operators in the Fock space. 
The Fock space can be defined by some orthonormal complete set of basis vectors: 

y^ \n){n\ = /, (n|r7i) = (5„„i. (3) 

n 

Let us introduce the additional Fock space which is identical copy of the initial Fock space 

(H| = 7, {n\m) = 5nrn- (4) 



^|n)(? 



We denote all vectors and operators in this additional Fock space by "tilde" . The vectors \n) and \n) span the so-called 
super- Fock space, which is a direct product of the original and the "tilde" Fock spaces. Operators in this super- Fock 
space will be called super-operators. Let us also introduce "left vacuum vector" (why it is appropriate to call it 
vacuum will be clear later) 



|/)=5^|n)®|n), (5) 



and " nonequilibrium wavefunction" : 

\pit)) = p{t)\I) =J2Prnnit)\n) ® |m), (6) 

nm 

where pmn{t) = {n\p(t)\m) . From Ti p{t) = 1 it follows that the left vacuum vector and the nonequilibrium wavefunc- 
tion are ortonormal to each other, {I\p(t)) — 1. With the use of the left vacuum and the nonequilibrium wavefunction 
we can re- write the average ^ as the following matrix element: 

{Ait)) = {I\A\pit)). (7) 

The Liouville equation ([T]) becomes equivalent to the time-dependent Schrodinger equation in the super-Fock space 

^l\pit))^L\p{t)), (8) 

where L ^ H — H is the Liouville superoperator (Liouvillian) . The left vacuum is always an eigenvector with zero 
eigenvalue of the Liouvillian (/| L = 0, which also automatically guarantees that {I\p{t)) = 1. 

Let us consider a system which consists of fermions distributed over N levels. Let us take vector \n) and \n) to be 

the particle number eigenstate \n) = |7ii?i2...njv) and \n) — |nlr^2■•■?^^f) 

alai\nin2...nN) = ni\nin2...nN)- (9) 

alai\nin2...nN) ^ ni\nin2...nN). (10) 

Here aj (a^) are fermion creation(annihilation) operators, which satisfy the standard anticommutation relations: 

{a^,a^j} = {ai,a'j} = Sij, {a^,aj} = {al,a^j} = {ai,aj} = {a-, a]} = 0, (11) 

and fermionic operators in initial and tilde Fock spaces anticommute with each other. Then the left vacuum can be 
written as 

|/) = 2J \nin2...nN) ® \nin2...nN) 

nin2---njv 

= exp{J24al)\0)^\0), (12) 

i 

where |0) and |0) are vacuums in ordinary and tilde Fock spaces (i.e. ai\0) = and ai\Q) = for all i). The vector 
\n) \n) can be always multiplied by some phase factor exp(iQ;). We found that, when one works with the fermionic 
systems, the convenient choice of the phase is the following 



l^)= E H)"^+"^+-+""l"i^2...n^)0|nin2...nw) 



nin2...nM 



= exp(-z^ala|)|0)®|0). (13) 

i 

One can readily demonstrate by the straightforward algebraic manipulations that 

a,\I) = -tal\I), al\I)^-ia,\I). (14) 

Following the terminology of the thermofield dynamics^^^ we will call the above relations as the "tilde conjugation 
rules" . The above relation which transform the original operators to the tilde operators is one of the most important 
relations and it is used very often in our derivations in this paper. In particular, it follows from (J14l) that the left 
vacuum in the super-Fock space, (/| , is the vacuum for a- ~ uij and a- + iaj operators. Since also (/| L = 0, it is 
appropriate to call (/| left vacuum vector. 

We will focus in this paper on the Lindblad quantum master equation.'*^ It can be written in the following general 
form 

|p(i) = -i[H, Pit)] + Y,i2L^pit)Ll - {it L^, Pit)}), (15) 



where H is the Hamiltonian of the system, L^ is a set of generally non-Hermitian Lindblad operators that represent 
the influence of the environment on the system, and {LtL^,/9(i)} means the anticommutator L^ L^p(t) + p{t)L'^ L^. 
The Lindblad master equation is the most general master equation which can be derived under the requirements that 
all probabilities are real and nonnegative, p{t) is always normalized, and p{t) can be obtained from p(0) by linear 
map.'**'^^ The last term in the Lindblad equation is so-called dissipator. It contains the term L^pL^ which creates 
quantum jumps between the states of the system and the term {LjjL^,p} balances the quantum fluctuations from 
the quantum jumps. 

Let us understand how we can re-write and solve the Lindblad master equation in super-fermion representation. 
We consider the following simple but nonetheless important example: A single level connected to thermal bath with 
temperature T and chemical potential p. The Hamiltonian is 

H = ea)a (16) 

and we take Lindblad operators in the following form^ 

ii = y/vla, U = x/f^al (17) 

The Lindblad equation becomes 

— p(t) — —ie{a*ap — pa* a) + Ti{2apa* — a* ap — pa} a) 
at 

-\-T2{'2,a* pa — aa* p — paa*) (18) 

If we act by this equation on the left vacuum |/), use the tilde conjugation rules (jl4[) and consider that the density 
matrix p = p{a\ a) is the operator in original Fock space therefore it commutes with all tilde operators, we obtain the 
time-dependent Schrodinger equation ([5]) for the nonequilibrium wavefunction \p{t)) with the Liouvillian L given by 

L = e{a^ a — a) a) — i{Vi — r2)(a^a + a^a) 

-2{TiZa + T25)a'^)~2iT2. (19) 

Within nonequilibrium thermofield dynamics a similar expression for Liouvillian (in nonequilibrium thermofield 
dynamics it is called Tildian) in super- Fock space was obtained axiomatically (see for example review^ ). 
Let us compute the average number of electrons on level e at time t: 

n(t) = {I\a^a\p(t)). (20) 

Differentiating this equation with respect to time we get 

z|n(t) = {I\a^aL\p{t)}. (21) 

Since (/| i = the right side of Eq. (PT|) can be written as the commutator {I\[a''a,L]\p{t)). Direct calculation of 
this commutator with the Liouvillian (TTOl) results into the following equation for the time-evolution of the occupation 
number 

^nit) = -2(ri + r2){n(t) - ^r^}. (22) 

at Fi + F2 

which has the solution 

n(t) = (no-noo)e-2(ri+r2)t + „^, (33) 

Here no is the initial occupation number a.t t — 0, rioo = p .^p — asymptotic occupation number at t = 00. Assuming 
IPoo) corresponds to the equilibrium density matrix in the grand canonical ensemble, we get 

(/|ata|p^) = /, (24) 

where f — [1 + exp[(e — p)/T]^^ is the Fermi-Dirac distribution function. So the choice of the rates Fi and F2 in the 
Lindblad operators in the form Fi = 7(1 — /), F2 = 7/ leads to asymptotic equilibrium state in the grand canonical 
ensemble with temperature T and chemical potential p. The parameter 7 = Fi -I- F2 is the relaxation rate to thermal 
equilibrium. The calculations along the same lines for the Lindblad master equation were also performed by Prosen 
within his "third quantization" methodi^ 

Let us summarize all our observation as a set of practical rules for the operations in the super-Fock space: 
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FIG. 1: Schematic illustration of the partitioning of the system. The left and right leads are assumed to be described by 
the non-interacting electrons, whereas electron-electron interaction is present in the molecule. Each lead is separated into two 
parts: The finite parts connected directly to the molecule and the macroscopically large parts (grey area) are represented by 
the Lindblad dissipators. 

1. The left vacuum vector |/) and the nonequilibrium wavefunction |p(t)) — p{t)\I) are invariant under the tilde 
conjugation \p{t)) = \p[t)), \I) = |/), and {I\p{t)) = 1. 

2. Tilde conjugation rules: The main rule is aj\I) — — ia^|/), aM) — ~ia\I) and as consequence the double tilde 
conjugation does not change the operator A = A and {ciA + C2B) — c^A + c^B. 

3. Evolution of the system is described by the time-dependent Schrodinger equation i-^\p{t)) — _L|p(f)), where the 
Liouvillian is obtained from the corresponding master equation for the density matrix with help of the tilde 
conjugation rules. The nonequilibrium average is given by {A{t)) — {I\A\p{t)) and {I\L = 0. 



III. QUANTUM MASTER EQUATIONS FOR ELECTRON TRANSPORT THROUGH THE 
INTERACTING REGION IN SUPER-FERMION REPRESENTATION 



Let us now apply the formalism of superoperators discussed in the previous section to the problem of electroir 
transport through the molecular junction. If we place a molecule into the contact with two electron reservoirs 
(usually metal leads) with different chemical potentials, the electric current starts to flow through it. Then, if wait 
for some time, which is much longer than typical electronic relaxation time of the molecule, the system will reach a 
nonequilibrium tiirre-independeirt steady state. In this section, we will discuss some general theoretical ideas of electroir 
transport through the molecule in nonequilibrium steady state regime based on the super-fermion representatioir of 
Lindblad master equation. 

We begin with the tunneling Hamiltonian: 



H ^ Hs{al,as)+ ^ ekaa\^ak 



k,a=L,R 



E 

s,k,a^L,R 



{tskaal^as + h.c), 



(25) 



where a 



it and 



ka i'^ka) creates (destroys) an electroir in the state k of either the left {a = L) or the right {a = R) lead, 
ai. ana as are creation annihilation operators in the molecule, and tgka are the coupling parameters between the 
system and the leads. In what follows we assume tgka — t*ska- Here Hs{al,as) represents the molecule and contains 
two-particle electron-electron correlations and, if necessary, electron-vibration coupling. The electron creation and 
annihilation operators aj (a^) and aj.^ (aka) satisfy standard anticommutation relations. We approximate each lead 
by N discrete single particle levels k coupled to the Lindblad dissipators (see Fig. [Ij . We assume that if two leads 
were disconnected from each other, they would be in thermodynamic equilibrium with temperatures Tl and Tn and 
chemical potentials fiL and fiR. As follows from the previous section, it can be accomplished by the use of the following 
set of 4N Lindblad operators: 



L 



kal 



V^kalO-ka, Lka2 — y'^ka2aka 



with Vkal = 7feQ(l - fka), Yka2 = Ikafka and fka = [1 + exp[(efca - tJ-a)/Ta] ^■ 

voltage changes only at the interface between the interacting region and the leads. 



(26) 

In the present model, the bias 
We also think that a particular 



choice of the dissipators is not an issue here as long as we can place the boundary between the part of the lead which 
is represented by the Lindblad dissipator and part of the lead which is given by the discrete single particle levels 
(Fig. [1]) deep enough inside the metal. This partitioning is justified in realistic systems since the screening length in 
the metallic leads is very short. 
The Lindblad master equation is 



dpjt) 
dt 

ka M=l:2 



= -l[H, p{t)] +Y.11 i^Lka^pmla^ - {^U^fc-A., PW})- (27) 



In principle, this Lindblad master equation (|27p for the electron transport can be obtained directly from the general 
full space tunneling Hamiltonian under some quite general assumptions by projecting out parts of the leadsi^^iiL 

Now we would like to convert this Lindblad master equation to the super-Fock space. Likewise to the single-level 
in equilibrium considered in the previous section, we act by the Lindblad equation on the left vacuum vector |J), 
employ tilde conjugation rules and use the fact that the density matrix p = p{a\a) is an operator in the original Fock 
space therefore it commutes with all tilde operators. Then the Lindblad master equation becomes the time-dependent 
Schrodinger equation in the super-Fock space ([5|) with the following non-Hcrmitian Liouvillian 



Here 



is the Liouvillian for the molecule and 



L = Lb+Ls + Lt. (28) 

Ls^Hs- Hs (29) 

Lt ^ -^^tska{a\^as-a\^s + h.c.) (30) 



is the Liouvillian, which describes the coupling between the molecule and the leads. They are both Hermitian. The 
Liouvillian for the leads 

is = ^ £ka{a\a°'ka " a\^aka) " i ^ LEfca (31) 

ka ka 

includes the non-Hermitian part of the Liouville operator which is responsible for the dissipation in the system: 

Hfca — (Ffccti — Tha2){aka^ka + O^ka^ka) — "^ii^kaiakaaka + ^ ka20.i.^aj^^) + 2Tka2- (32) 

It comprises the terms such as akaO,ka and a[.^a[.^ which usually appear in the theory of superfluidity.'^^ Due to these 
terms the structure of the nonequilibrium steady state wavefunction for the electron transport problem will have 
mathematical similarities with the vacuum in the Hartree-Fock-Bogoliubov theoryi^ This terms originate from the 
part of the Lindblad dissipator which is responsible for the "quantum jumps" between single-particle energy levels 
£ka in the leads. 

We are interested in nonequilibrium steady state situation, where the density matrix p[t) or nonequilibrium wave- 
function \p{t)) does not depend on time. Therefore, the electron transport problem is reduced to the problem of 
finding the eigenvector with zero eigenvalue of complex, non-Hermitian, finite-dimensional Liouville operator 

L\poo) = 0. (33) 

Here \poo) is nonequilibrium wavefunction for the nonequilibrium steady state. Subscript "cxd" in \poo) serves to 
emphasize that the nonequilibrium steady state can be also considered as asymptotic (t — > co) state of the system. 
Suppose wc have diagonalized the Liouvillian (|28p in some ^-modes exactly or within some approximation 

L = ^(S„e,U„ - K^iU- (34) 

n 

Suppose also that (/| ^j'j = (/| ^,\ = 0. Since {I\ L = 0, then after diagonalization L does not contain "c"-number terms. 
As a result a nonequilibrium steady state can be determined as a vacuum of ^- modes. The creation and annihilation 
operators of ^-modes will not be Hermitian adjoint to each other, since the left and right vacuums are different. If we 
recall that the basic idea of standard (equilibrium) quasiparticles is to represent true ground state of interacting many 
particle systems as a vacuum with respect to some quasiparticle annihilation operators ,^^'^" then these ^-modes can 



be regarded as nonequilibrium quasiparticles. The nonequilibrium quasiparticles include nonequilibrium, correlations 
and dissipation into their structure. This nonequilibrium quasiparticle description is, in principle, exact if one defines 
the quasiparticles in terms of the exact eigenstates of the many-particle Liouvillian. Although in practice it can 
be done approximately by establishing the relation between quasiparticles and "bare" particles of the system, for 
example, by canonical transformations. Below, we demonstrate how it can be done. 

Let us begin diagonalization of the Liouvillian ([28|) . First, we separate the part of the Liouvillian which contains 
only operators from the leads Lb- It is not Hermitian due to the presence of the dissipators. Direct calculations show 
that the following canonical (but not unitary) transformation 

Ofca — bka — ifkabf.^, ttka = bka + ifkabj^^, 

oL = (1 - /fca)^L + *^fe", 4a = (1 - /fea)^L ~ ^^ka, (35) 

with fka = ^ka2/(rkai + ^ka2) = [1 + cxp[(£fca - ^q)/Tq]"\ brings Lb to the diagonal form 

Lb=Y. {EkAjkc - E^JlJkc} . (36) 

ka 

Here Eka = £fca — i{^kai +^ka2) = £fca — ijka and E^^ = {Eka)* is complex conjugated energy. Nonunitary canonical 
transformations mean that these transformations, although being nonunitary, preserve the anticommutation relation 
between the fermion creation/annihilation operators. Since, b\^ = aj.^ — icika, the vector (J| is automatically the 
vacuum for 6[,^ and 6j,^ operators. We would like to repeat the important remark on the notation we use: only 
creation/annihilation operators written with letter a (such as, for example, Uka and a^^) are related to each other 
by Hermitian conjugation; all other creation (6^, c^, etc.) and annihilation operators {b, c, etc.) are "canonically 
conjugated" to each other , i.e. for example b'^ does not mean (6)t although {6,6^} = 1. 

Having reduced the reservoir part to the ideal gas of non-Hermitian fermions, let us begin to work with Liouvillian 
of the interacting region. We take the Hamiltonian for the molecule in the following general form 

f^s = Y1 Ksis2ali<^s2 + 4 H Vs,s2s:,s,al^al^as^as3, (37) 

SiS2 S1S2S3S4 

where ifsiS2 is single-particle matrix element which contains kinetic energy, electron- nuclei attraction and interaction 
with external fields such as, for example, gate voltage, Vs-^S2S3S4, is antisymmetrized matrix elements of Coulomb 
electron-electron interactions. Using the Wick theorem we perform the normal ordering of the Liouvillian for the 
interacting region: 

Ls=Hs-Hs = ii°^ + L's, (38) 

where L^ is the quadratic part 

^T = '^{I^s^s2 + Yl '^^4S3"^siS3S2S4) : a^^Os^ : -t.c, (39) 

S1S2 S3S4 

and Lg contains part of the electron-electron interaction, which is irreducible to the quadratic form : 

■^'s "" 4 Yl (^^I'ss^s^i : asi42«<*4a33 : -t-c) . (40) 

S1S2S3S4 

The steady state single-particle density matrix is ns2Si = {I\0'l-iO.S2 IPoo), and the notation (t.c.) means the tilde 
conjugation (i.e. a^ — >■ Og, rig^si ~^ "-s s j etc.). The normal ordering is asymmetric: it is performed with respect to 

the left vacuum (/| from the left and nonequilibrium vacuum |poo ), which is the solution of the following eigenvalue 
problem 

L^"^\P^^^)=0, (41) 

from the right. Here 

i(0) = L^°'> +Lt + Lb (42) 



8 

is the quadratic part of the Liouvilhan. We would hke to comment on the applicabiUty of the Wick theorem to our 
case, when the left vacuum state is not the same as the right vacuum. As it will be demonstrated below, we can 
always write operators as,al as a linear combinations of some fermionic creation and annihilation operators (let us 

cah them Cn ) in such a way that c„ |poo ) = Cn \poo) = 0, and (/|cjj = (/|cj, = 0. In this case, the Wick theorem is 
applicable^ and the Liouvilhan can be brought to the normal form with respect to the "left vacuum" (J| from the 
left and the vacuum for the quadratic part of the Liouvilhan |poo ) ((4T|) from the right. 

Since Lg is quadratic and Hermitian, it can be diagonalized exactly by unitary transformation D of creation and 
annihilation operators 



EDs 



'as' (43) 



Note, that D does not mix nontildc and tilde operators. To diagonalize pop . Lg — ^^. Esialag — aja^), the matrix 
D should satisfy the following equation 

J2{Ks^S2 + Yl '^^4S3^S1S3S2S4)^32S = ^S-DsiS- (44) 

S2 S3S4 

This is the nonlinear eigenproblem, since Us^sa depends on D and Ss- Next, we introduce system operators 
bl = al — icis, bl = al + ias which annihilate left vacuum (J| due to tilde conjugation rules. The canonically conju- 
gate annihilation operators arc bs ^ as, bs — as- Then the Liouvilhan takes the form 

L^Y. ^^^^^sbs -b%) + Lb + Lt + L's = L(0) + L'g, (45) 

S 

where the tunnehng interaction is expressed in terms of new creation and annihilation operators bka^ bs 

Lt = -J2 Tskc { [bljs + b%o.) - ibiJs + bfbka) + ifkc [blA + blJl)} , (46) 

ska 

and Tska — X^s' ts'kaDss' are renormalized matrix elements of the tunneling interaction. After the transformation 
the tunneling part of the Liouvilhan depends on temperatures and chemical potentials of the leads through the 
Fermi-Dirac occupation numbers fka- 

Next we diagonalize the quadratic part of p5|) . L^^\ exactly, 

L(") = ^(fi„4c„ - r!;2+H„). (47) 

n 

It results to the nonequilibrium Hartree-Fock theory for the electron transport problem. Creation and annihilation 
operators cjj, c„ can be regarded as nonequilibrium quasiparticles with complex spectrum J7„ (the spectrum of tilde 
conjugated nonequilibrium quasiparticles is given by — fi*). These operators obey the fermionic anticommutation 
relations, although c\ and c„ (c], and c„) are not Hermitian conjugated to each other. 

To find the internal structure and energy spectrum of nonequilibrium quasiparticles we use the equation-of-motion 
method; If the Liouvilhan is diagonal (|47p . then creation and annihilation operators c„ and c)^ must satisfy the 
following equations of motion 

:i;,LW]=-a„4, (48) 



c, 



[c„,LW]=a„c„. (49) 

Equations for c„, c]j are obtained from (|48I49|) by the tilde conjugation rule. We want to emphasize here that, since 
c„ ^ (cJi)^ and r2„ is complex, equations (|48I49|) can not be obtained from each other by the Hermitian conjugation. 
To avoid unnecessary complication, we do not give a general solution of Eqs. ipSl H^ . we would rather solve these 
equations for some particular examples in the next section. Here, we only note that operators which diagonalize the 
quadratic part of P5|) have the following form 



i = Y '^"' ^^^ + Y "^"^ '="^fc"' 

s ka 

n = ^(V'n^s&s +iipn,sb\) + ^ (V'n.fca&fca +iVn,kab\^), (50) 



ka 



and c^ and c„ are obtained from (|50)) by the tilde conjugation rule. Note that since the left vacuum (/| is the vacuum for 
6^, 6^ operators, nonequilibrium quasiparticle creation operators cjj, cjj are linear combinations of creation operators 
only. With help of the anticomniutation relations between nonequilibrium quasiparticle creation and annihilation 
operators, {cn,Cn,} = (5„„', {cn,Cn'} = 0, we find that the amplitudes tp, (p satisfy the following orthogonality 
conditions: 

/ ^ Yn, syn' ^ s ~r / ^ ^n, kaYn', ka — ^nn' 7 
s ka 

^(.i>n, sH>*n-, s - fn, sC', s) + ^{''Pn, kafn', ka ~ fn, ka^Pn' , ka) = 0. (51) 

s ka 

By diagonalizing the Liouvillian we simultaneously find the nonequilibrium steady state L^^'\poo) = as a (right) 
vacuum for operators c„ and c„. Using the transformation inverse to (j50p we can express any operator in terms 
of nonequilibrium quasiparticle operators. Then all physical quantities are calculated as an expectation value with 
respect to (/| and |pco ) vacuum states (see Eq. ([7])) in the nonequilibrium Hartree-Fock approximation. For the 
steady state current from the lead a = L, R 

k ks 

we derive, after some algebra, that 

-^a^ = -'>'^Tska{lpn,kafn,s " i'n,s'P*i,ka) ^ -"^ini^Tskalpn.kafn^s, (53) 

ksn ksn 

where the last equality follows from 

[bs,bua\ = ^{^n,s'Pn,ka " Vn,kaV*n,s) = 0. (54) 

n 

For the steady state single-particle density matrix we have 

SI = E ^si4£),,s^Vn,siV'«,4- (55) 



'3231 

Si ,s'r, .n 



Let us summarize the main results of this section. We have shown that the quadratic part of the Liouvillian ((28)) 
can be diagonalized by three canonical transformations. The first transformation diagonalizes the reservoir part of the 
Liouvillian including dissipator, the second is performed on the electrons in the interacting region and diagonalizes 
the Hartree-Fock part of Ls, and the last transformation mixes operators from the interacting region and the leads 
and diagonalizes the entire quadratic part of the Liouvillian. Two of these three transformations are not unitary, but 
all of them are canonical, which means that the anticommutation relations between particle creation and annihilation 
operators are preserved. The remaining part L'g can be taken into account via standard perturbation theory or 
nonperturbatively by, for example, coupled cluster method^^ or configuration interaction theory.^^ In fact, any other 
method for the correlated electronic structure calculations can be extended to nonequilibrium within our approach, 
although special care should be taken because one has to work with the nonunitary representation of the creation and 
annihilation operators. 

IV. ANALYTICAL AND NUMERICAL EXAMPLES 

A. Transport through a single-level molecule 

To illustrate the theory we consider a single-level molecule connected to two leads held at different chemical poten- 
tials. The molecule is described by the Hamiltonian 

Hs = ea^fa. (56) 

In this section we obtain nonequilibrium steady state wavefunction \poo) for this model and compare the results with 
the Landauer theory. To simplify the notation we assume throughout the calculations that the matrix elements of 
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the tunneling interaction are real and do not depend on the leads energy levels, i.e., tka — t^a ~ ^- ^*^ ^^^ Liouvillian 
becomes: 

L = e{a)a-o)a) - t'^{a\^a -a\^ + h.c.) +'^eka{a\^aka -a\oflka) -i^Hfea, (57) 

ka ka ka 

where Ilka is the dissipator for the lead a ~ L,R taken in the standard form (j32p . Our goal now is the exact 
diagonalization of the Liouvillian ()57p in terms of nonequilibrium, non-Hermitian quasiparticles: 

2N+1 

L^Y. {^ncW-K4.Cn)- (58) 

n=l 

To do this, we first perform the transformations (|35p over the leads operators and b^ ~ a^ — ia, b = a for the molecular 
level. Then the equations of motion ()48p give the following structure for the nonequilibrium quasiparticle creation 
and annihilation operators: 

ci=4'nb^+J2^n,kabl„, (59) 

ka 

Cn = -ipnb + ifnb^ + ^(V'n.fea^fca + i^n,kab\^). (60) 

ka 

The amplitudes ipm ipn,ka and quasiparticle energies fin are the solution of the following eigenvalue problems 

Si'n - t ^ i>n,ka = ^n'4'n, 
ka 

Ekai^n,ka " t^Jn = ^n'4'n,ka, (61) 

whereas the amplitudes (/?„ and ipn.ka satisfy the following nonhomogeneous system of linear equations: 

(e — ^,i)^Pn — ^2^ V'n.fca = ^Z_^ fka'>Pn,ka, 
ka ka 

{El^ - rtn)'Pn.ka - tlfin = -tfkalpn- (62) 

The amplitudes ipn, 4'n,ka should be normalized according to the first equation in ([5T|) . Operators c„, cj^ are obtained 
from ([59| and ([60]) by the tilde conjugation. Again, the transformations (|59I60I) are canonical but nonunitary, so 
{cn,c„/} = 5nn' and (c„)' ^ cjj, and nontilde and tilde operators anticommute. Although the analytical expressions 
for quasiparticle amplitudes and spectrum can be obtained (see Appendix), we found out that it is more convenient 
for practical calculations and for application to more complex system to solve Eqs. (pTl W^ numerically. Namely, we 
first solve the eigenvalue problem (1611) . and then with the known quasiparticle spectrum i7„ and amplitudes ^„, ipn,ka 
we solve the linear system of equations (|62|) . 

Since the Liouvillian (|57p is diagonal in terms of cj,, c„ and their tilde conjugate a^, Cn, we can associate the 
vacuum of c„, c„ with the nonequilibrium steady state: 

Cn\Poo) = Cn\Poo) = 0. (63) 

By the construction of c|i, cJ,, (/|c|( = (/|cj^ = 0. 

Using the transformation inverse to (|59[I60I) (see Appendix) we calculate the steady state current and nonequilibrium 
electron populations of the molecule and leads levels. For the current, according to ([53|, we get 

Ja = -2ilm^ V„,fca¥'n: (64) 

k n 

while for the population of the molecule and leads levels we derive 

{a) a) = (/| o)a \pao) = ^ V'n'iSn, 

n 
{a\^aka) = {I\a\^aka \Poo) = fka +2j^^.^a^n,fca- (65) 
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FIG. 2: (Color) Left panel: The current through the one-level molecule as a function of the energy of the level. We choose 
chemical potentials of the leads to be I-il,r ~ ±0.5 and Tl,r = 0.1. Right panel: The difference between calculated and exact 
currents for different values oi N = 100 x 2*° . 




FIG. 3: (Color) Left panel: The occupation number [Eq. (|65|l ] for the one-level molecule as a function of energy of the level. 
Right panel: The difference between calculated and exact occupation numbers for different values of A'^. The parameters are 
the same as in Fig. [2] 



As is expected, these occupation numbers are real (see Appendix). Moreover, using the explicit analytical expressions 
for amplitudes ipn,ka, 'fin,ka it can be demonstrated that the total number of electrons in both the leads is conserved 
quantity, i.e.. 



E' 



al^aka) 



E/^ 



(66) 



In what follows we consider the case of identical left and right leads, i.e., EkL — EkR — Ek- Then N eigenstates 
of (ICT]) coincide with the energies of the leads' levels, i.e., 51„ = Ei. Using analytical expressions for the amplitudes 
Tpn^ka and Lpn (sec Appendix) it can be easily shown that only eigenstates with J7„ = Ei contribute the sum over n 
in (j64p . (This result is also valid when i^ = atL-) It is similar to Meir and Wingreen observation that the expression 
for the current can be simplified within NEGF formalism when the spectral function of the coupling to the left lead is 
proportional to the spectral function of the couplings to the right leadj^ On the contrary, eigenstates with i7„ = Ei 
do not contribute to the molecule population in (I65p . 

Let us compare the results obtained within our approach with the results given by the Landauer theoryi^ The 
Landauer theory is exact for noninteracting electrons and it gives the following expression for the current^ 



1 

2^ 



duj- 






(67) 



and for the population 
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2ttJ {u-e-2A{Ljyp + {r{oj)y 



Here A(a;) and r(a;) are real and imaginary parts of the self-energy of the leads 

== A{uj) - - 



^H=t-E ,„_.i^,-. -AH^^rM. (69) 

fc 



Notice that the self-energy (|69p exactly appears in secular equation (IA3p . We do not distinguish between the self- 
energies of the left and right leads, since the both leads are identical. In numerical calculations by the Landauer 
theory we put A'^ — >■ oo, 7 — > and assume that leads densities of states are constant within the energy bandwidth 
[Emin '■ Emax\- Then T{uj) = r for Efnin < '^ < Emax and zero otherwise. For K{uj) we have the Cauchy principal 
value integral 



CXD 



T{e)de _ T 

id — e 2tt 



uj — E„ 



(70) 



—00 



In our calculations with the kinetic equation we assume that N energy levels in each lead are evenly spaced in the 
bandwidth [Emin ■ Emax] = [—5 : 5]. The tunneling coupling strength t is computed from the F = 2TTrit'^ — 1, where 
r] = N/(Ejnax — Ejnin) IS the density of states. Below we put 7 = 2Ae, where Ae is the energy spacing between states 
in the leads. Under this choice of the parameters F(a;), as is computed by (p^ . is constant equal to 1 within the 
bandwidth of the leads and vanishes when uj is outside the bandwidth. 

In the left panel of Fig. [21 we plot the current [Eq. (|64l) ] calculated within our approach for N — 100 as a function 
of the level energy along with the exact current computed by the Landauer formula (j67| . We see that the calculated 
current agrees well with the exact one. The largest deviation between the two currents is obtained when the currents 
reach their maxima. In the right panel of Fig. [2] we plot the difference between the calculated within our approach and 
exact currents for different values of N. It is evident from the figure that the difference becomes smaller as the leads 
densities of states increase. It agrees with our observation (see Appendix) that our expression for the current ([M)) 
becomes the Landauer formula (|67p in the limit of macroscopically large leads. We also calculate nonequilibrium 
population of the molecule and compare them with the exact values (p5)) for different values of A^. The results 
are shown in Fig. [3l Again, we reproduce the exact results by increasing the part of the leads included into the 
consideration. 

B. Out of equilibrium Anderson model in Hartree-Fock approximation 

In this section, we apply our method to a physically more interesting example, namely to electron transport through 
a spin-degenerate single level with local Coulomb interaction (so called Anderson model). The Hamiltonian for the 
molecule has the form 

Hs ^ e'^alaa + Ua'^^a'^alai, (71) 

where aj., a^. are the creation and annihilation operators for the spin- up {a =t) and spin-down {a =4,) electrons in 
the molecule, e is the energy of the single level in the molecule, and the charging energy U characterizes the Coulomb 
interaction of electrons in the molecule. Here we restrict our consideration by a nonmagnetic system, therefore the 
energy levels in the leads are spin degenerates and all the coupling strengths tkaa- take the same value t. Under these 
assumptions we have 

Hb+ Ht = '^ekaa\^^akaa - t'^{a\^a°^c, + h.c.). (72) 

kaa kaa 

We are interested in quantum transport in Coulomb blockade regime, where we charging energy U is much greater 
than the effective coupling F^^^ = 2'K'qt^ between the molecule and the leads. So we work in the regime of strong 
Coulomb interaction and weak coupling to the leads. 
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FIG. 4: (Color) The current (left panel) and the occupation numbers (right panel) for out of equilibrium Anderson model at 
fixed gate voltage Vg — —0.2 as a function of voltage difference between left and right leads. For the current we show the 
results obtained with different number of states A'^ in each lead and the exact Hartree-Fock result. In the right panel the total 
occupation number and the nonequilibrium populations of the Hartree-Fock levels are shown for A'^ = 2000. 



Using the Wick theorem with respect to nonequilibrium steady state vacuum \poo) (which is yet to be defined) 
from the right and the left vacuum (/|, we obtain the Hartree-Fock part of Hg in the following form 



H% = ^(e + Un^„)a\a„ - [/(n^^aja^ + ntiajaj,) 



(73) 



where Ua = (/| a\aa |poo ) and n„„i = {I\ a^^,a„ \Poo)- Since the system is nonmagnetic, the occupation numbers for 
both spin orientations are the same n„ — n^„. Moreover, due to the symmetry reasons n^^ = n^^. Then the following 
unitary transformation 



vt 



;^(4-«l) 



„t- 



-^(4 + 4) 



brings i?S to a diagonal form 



2 - ^V-t ' "V 



(74) 



(75) 



(76) 



where the Hartree-Fock mean-field energies are 

£1(2) = e + U{7ia ± n^i) = e + Un2{i), 

and Hg = (/| oi\as \poo) is the occupation number for the Hartree-Fock single-particle level "s" (the crossover matrix 
element (/| a\as'^s \poo) is equal zero) and s — 1(2) corresponds to " + (— )" signs respectively. We can regard a\ as 
the antibonding orbital creation operator, and a^ as the creation operator for the bonding orbital in the dot. From 
the above equation it follows that if n^^ ^ 0, the Hartree-Fock bonding/antibonding energies £1(2) are split and the 
populations of bonding/antibonding orbitals are generally not the same ni 7^ 712. 

Following the procedure described in Sect. Hill we introduce h\ and h\^^ operators which annihilate the left vacuum 
(J| . Then the quadratic part of the Liouvillian becomes 

L(")=^£,(fet6,-fet5^) + LB 

S 



(77) 



kaas 



where Ti 



-t/^/2 and all other T^^. are equal t/\/2. Likewise the noninteracting single-level model considered 
in the previous section, the Liouvillian (j77p can be diagonalized in terms of nonequilibrium quasiparticles by the 



U 
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FIG. 5: (Color) The current (left panel) and the occupation number (right panel) for out of equilibrium Anderson model at 
small applied voltage fiL,R ~ ±0.025 as a function of gate voltage. The parameters and notations are the same as in Fig. ID 



equation of motion method. The structure of nonequilibrium quasiparticles, which diagonalize (|77|) is the following 
(n=l,...,2 + 4iV) 



s kcxa 

s kaa- 

where the aniphtudes ip^ ip and the spectrum Q^ cire obey of the equations (s = 1, 2) 

kaa 



(78) 



(79) 



and 



kaa ka 

\^ka ~ ^ ^7i)^7i,kaa ~ / ^ ^ sa^n.s — ~Jka / ^ '^saWn 



(80) 



The amphtudes depend on Hartree-Fock occupations numbers ris, which in turn depend on the amplitudes: 






V'n.s'Pn, 



(81) 



Therefore the equations (|79I80I) should be solved via self-consistent iterations. After the diagonalization, the Liouvil- 
han (j77p takes the form (j47p . Having found the amplitudes for the nonequilibrium quasiparticles we can compute the 
nonequilibrium population of the molecule {n-^ + rt^ = ni + 112) as well as the steady state current ([53|). 

Likewise the example with single molecule level, in case of identical leads the system (j79|) have N twice degenerate 
eigenstates equal to the energy levels in the leads. Only these eigenstates contribute to the steady state current, but 
they do not contribute to the population of the molecule. 

In our further numerical calculations we assume the symmetrical voltage drop, /xl,_r = Ep ± 2^' ^"^'^ ^p ~ '^• 
Also, we introduce gate voltage Vg, so that the energy of molecular level is e{Vg) = Ep + Vg. We put U = 1.0 for 
the strength of the Coulomb interaction. The leads energy levels are spaced in the bandwidth [-1:1] and r-|- 4^ = 0.01, 
Tl,r = 0.001. 

Fig. SI shows the calculated current through the molecule and the occupation numbers as function of the bias 
voltage applied across the junction at fixed gate voltage Vg = —0.2. While both chemical potentials are above the 
molecular energy level e = —0.2 the electric current through the molecule is small and the average occupation number 
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of the molecule is close to unity. It means that the bonding Hartree-Fock level s = 2 is below ^jiL,R {^2 = s) and it is 
totally occupied (n2 = 1) while the antibonding level s = 1 is above piL,R (£i = e + U) and it is empty (ni = 0). The 
first step in the current occurs when the right chemical potential reaches the energy e {V ~ 0.4) and the Hartree-Fock 
bonding level s = 2 becomes involved in electronic transport. As a result the occupation of this level decreases to 
0.5, and the energy of the empty Hartree-Fock level s — 1 reduces to ei = e -I- 0.5C/. The total number of electrons is 
the molecule becomes ni + n2 — 0.5. The second step in the current occurs when the left chemical potential reaches 
the energy level ei = e -I- 0.5[/ {V ~ 0.6) for the antibonding orbital and electrons starts to flow through it. This 
increases the occupation number of the antibonding Hartree-Fock orbital level s = 1 from zero to 0.5 and, therefore, 
the bonding level s = 2 is pushed up from £2 = £ to £2 = e + Q.bU. At higher voltage {V > 0.6) the Hartree-Fock levels 
become degenerate and they both lay between /i^ and /i^. Since each of them is half-occupied, the total number of 
electrons in the molecule becomes again equal to unity. 

Fig. [5] shows the current and the occupation numbers as a function of the gate voltage Vg, at small applied voltage 
IJLl,r = ±0.025. At Vg < —1.0 both Hartree-Fock levels are below (ei_2 — e + U < 0) the chemical potentials for both 
leads. Accordingly, the occupation number is equal to two and the current through the junction is close to zero. At 
Vg « —1.0 the molecule starts to conduct electrons This result in a pronounced increase in the current accompanied by 
a decrease in the occupation number. At — 1.0 < V^ < the Hartree-Fock bonding level is below chemical potentials 
(£2 — e), and antibonding level is above them (ei = s + U). Therefore the occupation number is equal to unity and the 
current through the molecule becomes almost zero At Vg « 0.0 the bonding level becomes involved in the electronic 
transport giving rise to the current increase. At higher gate voltage, when the single-particle level e is above the 
chemical potentials, the molecule become empty and the current again drops to zero. 

In Figs. |3]and[S]we compare our results with the exact Landauer current through Hartree-Fock levels. To obtain 
the exact current we first solve iteratively the nonlinear equations ?^l(2) = ?^ea:act (£1(2)) (|68p . Then the Hartree-Fock 
levels £1(2) are used in (j67p to compute the current. Figs. U and [S] show that the currents calculated within the present 
approach converge to the exact result with increasing value of N . 

V. CONCLUSIONS 

Based on the super-fermion representation of the Liouville space, we developed an approach which enables us to 
transform the quantum master equations to the super-Fock space and then to use standard methods of quantum 
field theory to solve them in the nonequilibrium steady state regime. We worked with the Lindblad master equation 
in this paper although the derivations can be readily extended to more sophisticated master equations. The main 
technical difficulty of the approach is that the left vacuum is always different from the right vacuum state unless we 
have a system in thermodynamic equilibrium. This prohibits the use of the unitary transformations to diagonalize 
the Liouvillian. The problem was circumvented by the development of a set of nonunitary, canonical transforma- 
tions between particle creation and annihilation operators. These nonunitary canonical transformations preserve the 
anticommutation relations between the fermionic operators and significantly facilitate the derivations. 

Starting with Lindblad master equation for the electron transport through the interacting region, we converted the 
problem of finding the nonequilibrium steady state to the many-body problem with non-Hermitian " Hamiltonian" in 
super-Fock space. Then we demonstrated that despite the fact that the left vacuum is different from the right vacuum 
we still can use the Wick theorem. Using the Wick theorem we transformed the Liouvillian to the normal ordered 
form, introduced nonequilibrium quasiparticles and developed a general many-body theory for electron transport 
through interacting region. We applied the approach to electron transport through a single level molecule. Then we 
considered the system with electron-electron interactions, namely out of equilibrium Anderson model in Hartree-Fock 
approximation. The Wick theorem was applied to obtain the Hartree-Fock solution of the transport problem. We 
demonstrated that it is consistent with the Landauer theory, which is exact for these models. Being formulated in 
the language of Fock spaces, creation and annihilation operators and normal ordered "Hamiltonian" the proposed 
approach is not only capable of doing perturbative calculations but also has a great complementarity to the nonper- 
turbative many-body methods of molecular electronic structure calculations such as, for example, coupled cluster or 
configuration interaction theories. 
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Appendix A: Structure of nonequilibrium quasiparticles creation and annihilation operators for single level 

model 

For the single level model the analytical solution of Eqs. (|61] l62| has the following form: 

V'n=AC'/', ^n,kc. = -^ ppV-n, (Al) 

rLhn — iir. 



Kd ^^n 



•s—\ Jka^ka 



Vn = A/„ ' '^ , ipn,ka = -^^ T^K'Pn ^ l/jnfka) ■ (A2) 

7feQ ^ka ^ ^'n 



ka i^ka - ^n){El^ ~ fin) 

where Mn = 1 + t^ X](-^fca ~ ^n)^^, and fi„ obey the secular equation 

kct 



kOL 

Some constraints on quasiparticle amplitudes can be obtained by using the transformations inverse to (|59l60p : 

n n 

^L = X! '^".fc"'^"' ^fc" = ^{'^n^kaCn - i^l^kcA) (A4) 

n n 

(the expressions for h, h\ b^cn b\^ are obtained from (|A4[) by the tilde conjugation). In particular, from {6^0,, bka] — 
and {6, &} = it follows that ^ '4'n,kafn,ka and ^ V'n'ySTi are real. 

n n 

If we have identical left and right leads, i.e., if E^l — E^r — Ek, then N eigenstates of (|6T|) coincide with the 
energies of the leads levels, i.e., ri„ = £';. For such normalized eigenstates we have 

V'« = 0, 1pn,kL^-'4'n,kR^Skl—l=, (A5) 

v2 
where Ski is the Kroneker symbol, and 

^«-2-V. ^(f^^-f^-^ , ^„,. = -^. (A6) 

fe i-fc'i -%■) 

This result is also valid when Ir = at^. The only difference is that ipn.iR = —ct^'njL- Therefore, in case of identical 
leads the summation over n in (|64p can be divided into two parts. The first part involves summation over A^ + 1 
solutions of Eq. (jA3l) . while the second part involves summation over N solutions of (|5T|) such that rj„ = £";. The 
first part does not depend on index a and it vanishes. Therefore 

J ^j ST ^^(/jL - fm) ,.„. 

J a = Tim 2^ -^ . (A7) 

Here the upper (lower) sign corresponds to a = L [a = R). It is obvious that in the limit of macroscopically large 
leads [N — > 00, 7—^0) eq. (jA7p becomes the standard Landauer formula (l67l) for the current. 
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